Reconstructing large running-index inflaton potentials 
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Recent fits of cosmological parameters by the first year Wilkinson Microwave Anisotropy Probe 
(WMAP) measurement seem to favor a primordial scalar spectrum with a large varying index from 
blue to red. We use the infiationary fiow equations to reconstruct large running-index infiaton 
potentials and comment on current status on the inflationary fiow. We find previous negligence 
of higher order slow rolling contributions when using the flow equations would lead to unprecise 
results. 

PACS number(s): 98.80. Cq, ll.lO.Kk 



In the past decade, inflation theory has successfully 
passed several nontrivial tests. In particular, the re- 
cently released Wilkinson Microwave Anisotropy Probe 
(WMAP) data Jj have detected a large- angle anti- 
correlation in the temperature-polarization cross-power 
spectrum, which is the signature of adiabatic superhori- 
zon fluctuations at the time of decoupling^. The high- 
precision WMAP data have been used to fit the cosmo- 
logical parameters and confront the predictions of in- 
flationary scenarios respectively in Refs. |3 and 0. It is 
noted that there might be possible discrepancies between 
predictions and observations on the largest and smallest 
scales. Basing on this fact Spergel et al introduced 
a new parameter: dnsjdhik in the fit. With the global 
fit to WMAPext, 2df H and Lyman-a forests "s^the au- 
thors found a preference of nonzero running at around 
2(7. Although the use of Lyman-a data was questioned 
by 0, the fit to WMAP alone favored the running to 
about 1.5a: at the pivot scale ko — 0.002 Mpc~^, the 
best-fit values of the scalar power spectrum for WMAP 
alone being = 1.20^0 and dus/dlnk = -0.077i|^:g^^ 
13 . It is recently found that the inclusion of VSA data 
can lead to a negative running at a level of more than 
95% confidence when fitting to CMB only (3. Theoreti- 
cally there have been studies in the literature since the 
release of the WMAP data on models of inflation which 
provide a running index required by the WMAP |^ . 
On the other hand, Bridle et al claimed that the need 
of running should come from the low CMB quadrupole 
alonejlOj and thereafter many theoretical models were 
put forward Meanwhile Mukherjee and Wang 
used the model independent reconstruction of primordial 
spectrum and found that a running index from blue to 
red was indeed favored. 

The need of large running has been studied widely|^ 
ES in III 13 after the release of flrst year WMAP data. 
At least WMAP could not rule out significant running 
and if further stands, it could severely constrain inflation 
model buildings. In the present paper we will use the 
inflationary flow equations to generate potentials which 



can lead to large variations in the spectral index. 

The inflationary flow equations were firstly introduced 
by Hoffman and Turner |l5j to study generic predictions 
of slow-rolling inflation and for the first time compared 
with observations in RefEl Latter Kinney 0| developed 
the equations as the inflaton potential "generator" , which 
has been used widely in the literature 0, Q, It was 
later shown by Ref.[l9j that the flow equations could be 
used to reconstruct inflaton potentials. By randomly gen- 
erating a large number of initial flow parameters, evolv- 
ing towards the observed scale and comparing with the 
observations, inflaton potentials that satisfy the obser- 
vations are then reconstructed |19|. Thereafter Liddle|20j 
pointed out subtly some disadvantages with current ver- 
sion of inflationary flow. We will reconstruct inflaton 
potentials with large running of spectral index, comment 
on Liddle's remarks and put forward our understandings 
on inflationary flow. 

We first follow closely the notations by KinnevfT^. For 
a single scalar field inflaton (p it obeys the equation of 
motion 



4> + 3H^ + ¥'{((,) = 0, 



(1) 



where H = (a/a) is the Hubble parameter and V{4>) is 
its potential. Another equation is the Hamilton- Jacobi 
equation 
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The scale factor during inflation is given by 
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where the number of e-folds N is 
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The slow-roll parameter e can be expressed by 
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in the convention of Hubble expansion. Accordingly, high 
order slow-roll parameters in Hubble expansion can be 
given |2ll |: 
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it is convenient to take the derivative with respect to 
number of e-folds instead of (/). The evolution of above 
slow-roll parameters can be described by the "flow" equa- 
tions |ilE3 
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As shown by R.eflT^ the flow equations have a stable late 
time attractor with 



e = 'Ah - 0, 
(T = const. 



(10) 



For any single field inflation model with common dy- 
namics (i.e., satisfying above equation of motion and 
Hamilton- Jacobi equation ), the background dynamics 
of the Hubble parameter (up to a normalization) can be 
equivalently described by flow equations as long as I can 
be extended to infinity. However, numerical calculation 
cannot accommodate infinite set of equations and one 
has to truncate above flow equations. A common trun- 
cation to M — th order assumes the ^Ah are all zero for 
I > M. After truncation the flow equations cannot ac- 
commodate all the single field inflation models, as shown 
by Liddlej2fl,], but the flow equation itself is still exact. 
We have 
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For given initial values of </> and H, their trajectories can 
be correspondingly given via the flow equations, V{4)) can 
then be worked out with the Hamilton- Jacobi equation 
once H is known. The initial value of 4> is arbitrary and 
one can set it to be zero, but H has to be determined on 



the observable scales from the primordial scalar spectrum 
(e.g., from WMAP ^ ) : 
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The inflationary "observables" include tensor/scalar 
ratio r, the spectral index rig, and the "running" of the 
spectral index dns/dlnk, etc. To the second order in 
slow roU(SR) one gets HlHi, 



r = 16e [l-C{cT + 2e)] 
for the tensor/scalar ratio, and 



(13) 



ns-l = (T-(5 - 3C) e^-i (3 - 5C) ae+i (3 - C) (^Ar) 

(14) 

for the spectral index. C = 4(ln2 -I- 7) — 5, where 7 ~ 
0.577 is Euler's constant. dns/d\nk can be given via the 
relation 
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which is tedious but can be directly worked out from 
Eq. H14|l . The number of e-folds N which corresponds 
to the observable scale is largely uncertain due to the 
uncertainty in the energy de nsity during inflation and 
the reheating temperature [2^ |24| . typically in the range 
[40,60]. 

Despite the doubt on whether flow equations corre- 
spond directly to inflationary dynamics, we can recon- 
struct inflaton potentials which lead to large variations in 
the spectral index using the Monte Carlo method, which 
may be still exact after a truncation at higher order in 
the flow equations. Similar to RefHol the algorithm for 
our Monte Carlo reconstruction is as follows: 

1. Specify a "window" of parameter space of the ob- 
servables: e.g. 1(7 WMAP constraints on — 1, r 
and dns / d\nk and their associated error bars. 

2. Select a random point in slow roll space, [e, (5, ^Ah], 
truncated at order M in the slow roll expansion. 

3. Evolve forward in time {dN < 0) until either (a) 
inflation ends (e > 1), or (b) the evolution reaches 
a latc-timc fixed point (e = ^Ah = Q, <y — const), 
or (c) I^AhI > 100, or (d) inflation docs not end 
after evolving AiV > 1000 . 

4. If the evolution reaches a late-time fixed point, goto 
6 (Our main intention is to search large dns/dlnk, 
which is zero is this case ) . If the evolution reaches 
I^AhI > 100, for the sake of saving computing time 
and to ensure the validity of expanding the observ- 
ables to second order in slow roll, goto 6. If inflation 
does not end after evolving AN > 1000, record it 
as insignificant point, goto 6. 
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5. If inflation ends, evolve the flow equations back- 
ward 40 until 55 e-folds from the end of inflation. 
Calculate the observable parameters at each point. 
Once the observable window is satisfied, record it 
as nontrivial point, compute the potential and add 
this model to the ensemble of "reconstructed" po- 
tentials and go to 6 . Else evolve until 55, calculate 
the observables and record as a nontrivial point. 

6. Repeat steps 2 to 5 until the desired number of 
models have been studied. 

We truncate the slow roll hierarchy to 5th order, the 
model parameters are randomly from the following uni- 
form distributions: 

e = [0,0.8] 

0- = [-0.5,0.5] 
^Ah = [-0.05,0.05] 
^Ah = [-0.025,0.025], 

'Ah = 0. (16) 

The thinning on the initial slow roll hierarchy is set to en- 
sure a reasonable convergence. However it is still possible 
that one gets larger values for higher order parameters, 
so we set a "lock" that |^Ah| < 100 in our simulation. As 
we are merely intending to find several potentials which 
can give large dns / d\nk, our method is in this sense ap- 
plicable. For a 1000,000 simulation, we get 

• Latc-time attractor: 921,793. 

• Nontrivial points: 75,300. 

• I^AhI > 100: 2,906. 

• Insignificant: 1. 

A companion running of truncation to 8th order is also 
tried for crosscheck. In this case we use almost the same 
method as the WMAP team 01 except the negligence of 
I^AhI > 10^ (this is taken as Method Crosscheck below). 
This time we set the range N = [40,70], we get (for a 
1000,000 simulation) 

• Late-time attractor: 924,164. 

• Nontrivial points: 75,646. 

• I^AhI > 10^: 187. 

• Insignificant: 3. 

We plot our zoo of the nontrivial points in Fig.l. In- 
terestingly there are some straight lines for n ~ 1.01 and 
dns/dlnk ^ —0.02, we shall make further comments on 
this below. Meanwhile the strait lines are unavailable 
for Method Crosscheck, as shown in Fig. 2. The "win- 
dow" we open here is around la WMAP constraints 
at A: = 0.002 Mpc'^: 1.01 < ns < 1.4, r < 1.14, 



—0.021 < dns/dlnk < —0.13. For the primordial spec- 
trum with large running in the spectral index, it is ex- 
tremely difficult to make an exact description on the re- 
sulted spectrum: WMAP team assumed a constant run- 
ning on spectral index in their fit and gave the above 
1(7 region. However from the evolution of flow equations 
dns/dlnk is generically not a constant. Under such cir- 
cumstances many wrong conclusions would be reached 
once only above window is used. Firstly small nonzero 
d{dns / dink) / dink other than the assumption taken by 
WMAP team is no doubt a possibly good fit to WMAP, 
the observational data can never restrict higher order 
terms to be exactly zero. Secondly what WMAP con- 
straints on is a continuous region from 10"** ^ 6 x 10~^ 
h.Mpc~^, the parameters from fiow may change signif- 
icantly within few number of e-folds, this may lead to 
some obvious mismatching. Theoretically flow equations 
can even be used to rebuild potentials which lead to sup- 
pressed large scale primordial spectrum or those which 
get a kick on smaller scales to achieve larger CMB TE 
multipoles on the largest scales. But the window is ex- 
tremely difficult to set. In the case of a large dns/dlnk 
we have to make some evaluations by hand. The first step 
is that we open the "window" as stated above, we get 858 
points which appear in the "window" for the 1000,000 it- 
erations and 8,788 points for a 10,000,000 simulation. We 
plot the 8,788 "raw" data in Fig. 3. We make a hand se- 
lection around the first 100 "raw" points, the data with 
global blue n and smaller running have been extracted. 

After careful extracting we find 81 of the data lead to 
global blue index for > 30 (when reanalyzing N is 
loosened to 30, which is still acceptable respecting the 
reheating temperature limit ), 5 where the constraint on 
ns and dns/dlnk cannot be simultaneously satisfied, 2 
where second order SR does not work well (for detailed 
discussions see below), only 12 points are left, leading 
to running enough from blue to red (It is theoretically 
plausible to exclude those which lead to global blue in- 
dex, which could be added to the window). We show 
the 12 resulting ns in Fig. 4, their initial flow parameters 
are shown in Table 1. The corresponding potentials are 
shown in Fig. 5 and their trajectories in Fig. 6. The tra- 
jectories are rather complicated. Very interestingly some 
trajectories overlap during some period in Fig. 6, as can 
also be seen from Fig. 5, around the CMB and Large Scale 
Structure interest scale (black/dark lines) some potentials 
do overlap, but are then separated on other scales. We 
find that r is no more than 0.5 for the twelve points. Sim- 
ilarly, if necessary, one can work on all the 8,788 "raw" 
data and filtering out hundreds of potentials that lead to 
large variations in the index. 

There have been in the literature several works which 
try to reconstruct inflaton potential using ideal observa- 
tional data^i lU El]' Generally speaking, flrst year 
WMAP data have not provided stringent window for 
one to constrain the inflaton potential. WMAP has 
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FIG. 1: Models plotted in the {ns,r) and {ns,dns / d\nk) 
plane for an M = 5 Monte Carlo. The cause of short straight 
lines around ng = 1 and dns/dlnk = —0.02 can be found in 
Fig.7. 



FIG. 3: Models plotted in the {ns, dns/dlnk) plane for 
a 10,000,000 simulation fitting into the window around la 
WMAP constraints at fc = 0.002 Mpc"\ The cause of the 
straight line around ns ~ 1 can be found in Fig.7. 




FIG. 2: Models plotted in the (ns,r) and {ns, dns/dlnk) 
plane for an M = 8 Monte Carlo. The number of e-folds is 
now randomly generated in range [40, 70]. 



not detected tensor contributions in the primordial spec- 
trum: r < 10. This leads to e < 0.07 and hence 
< 4.8 X 10-i°m|,;, this gives < V" < 5.8 x 10-"m|,; 
via the Hamilton-Jacobi equation. As also shown al- 
though a large running is favored, WMAP is consistent 
with scale invariant spectrum. Peiris et.al.0 categorized 
the slow roll parameters into four classes for generic single 
field inflaton potentials and found for either class there is 
broad region consistent with current observations. In this 



sense first year WMAP data (together with other obser- 
vational data) have not been able to work as a stringent 
enough discriminator. As for the reconstructing method, 
there would be broad parametric space when the obser- 
vational window is less stringent. However as the ob- 
servational data accumulate and the error bar shrinks, 
e.g. when tensor contribution is exactly measured and 
running of the spectral index is strongly confirmed, re- 
constructing the favored inflaton potential would make 
its sense. On the inflationary flow itself, it would be in- 
efficient since very small fraction fall into our window 
above, meanwhile most of the iterations fall into lat-time 
attractors and little part satisfies the window. In any case 
it provides a testable inflaton potential generator which 
may be exactly solvable. Liddle|20j has made some de- 
tailed descriptions on its shortcomings: Firstly the flow 
equations (Eq.9 in this paper) do not relate directly to 
the inflationary dynamics. This should not be a severe 
problem since the reconstruction of potential needs using 
the Hamilton-Jacobi equation, where the dynamics is in- 
cluded. Secondly due to the truncation on M — th order, 
the effective potential is only in such form 

y(0) = ^H^l + A^c^+.-. + AM+ir^'Y (17) 

X \ 1 "^ii Ml + ■ ■ ■ + {M + i)AM+i^^' y" 

^ 3 An \1 + Ai(P+--- + Am+i(I)^'^+^ J 

We made a detailed check and found the form is satis- 
fied exactly. This in fact does accomodate many inflaton 
potentials as there are so many undecided parameters to 
choose and an ideally no truncation would accomodate 



FIG. 5: Twelve reconstructed inflaton potentials from Fig. 4. 
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FIG. 4: Twelve resulting spectral indices that runs signifi- 
cantly from blue to red, satisfying la WMAP constraints^. 



all the single field inflaton potentials which satisfy above 
equation of motion and the Hamilton-Jacobi equation. 
A main loophole is pointed by Liddle that for the 
nontrivial points where inflation ends at e = 1 the full 
potential is negative at its minimum. For example Eq.l7 
reduces to V{(j)) = Vo{{l + Ai^)^ — when only e is 
nonzero where the flow equation starts to evolve (we set 
4> to be zero here). However as the flow equation only 
describes the dynamics of inflation at e < 1 where V is 
always positive (as can be clearly seen from Eq.2 ), we 
can assume that the exact potential after inflation can 
be matched by other potentials which reach the mini- 
mum at 1^ = (or arguably matched at e = 1), the 
perturbations would certainly be the same as in Eq.l7. 
In this sense flow equations can reconstruct potentials 
during inflation. Similar arguments hold for the fixed 
points where inflaton reaches a local positive minimum 
and drives eternal inflation, this situation can be related 
to the standard hybrid inflation mechanism l^^, as 
shown by Easther and Kinney case the flow 
equations seem to be inefficient as the inflaton potential 
generator when confronted with a large negative running 
from blue to red, which is around the la region of first 
year data of WMAP. Further investigation to find bet- 
ter infaton generator is necessary confronted with future 
precise observations. 
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FIG. 6: Twelve trajectories of e-a from Fig. 4 during inflation. 



So far we have left the second order SR approxima- 
tion(SRA) undiscussed (We thank the anonymous refer- 
ees for inspirations on this issue) . It has been noted in the 
literature 28J that higher order SR contributions may not 
be negligible especially when confronted with high preci- 
sion observational data. This being the first time to check 
the precision of SR approximation in the framework of 
inflationary flow, we first rewrite the spectral index to 
the third order in the form of flow parameters (Basing on 
RefsIHIlS): 

ns-1 = cr + 2a^AH + ( a'^ - ^Ah + (4 - 12a) 
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FIG. 7; Crosscheck of SR approximations for a 10,000,000 
simulation on the observable values of ns and dns /dink, with 
dns jd In k being the vertical axis. The left panel shows those 
which lead to wrong observable values due to nonegligible 
higher order effects of 'Ah (' > 3); The middle and right panel 
show those which give wrong values of ns and dns /din k due 
to insufficient SR approximation to second order. 
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FIG. 8: The same as Fig.7 for a 10,000,000 iteration on 
Method Crosscheck described in the text. 
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where a = ^^-^ = 0.7296. We made three sets of cross- 
check in our 1000,000 and 10,000,000 iterations. Firstly 
one can see from Eq.l8 the factors on ^Ah and '^Ah differ 
no more than five times. There being no efforts avail- 
able in higher order SRA as Reffi^ has explicitly did, 
we first make a strong limit that maxi^ \ '^Ah |, | ^Ah |, 
I ^Ah I) < 0.01 (taken as Condition A)is satisfied. Those 
which violate this condition are plotted in the left panel 
of Fig.7. Other limits are taken to ensure that the val- 
ues of dns /dink (taken as Condition B)and ns (taken 
as Condition C)differ no more than 0.01 when a third or- 
der SRA is assumed instead of second order SRA. They 
are delineated in the right and middle panels of Fig.7. 
To our great surprise 319, 134 and 6376 points violate 
Conditions A,B and C respectively in the 1000,000 it- 
eration and for the 10,000,000 iteration 3089, 1404 and 
64589 points violate the three conditions ( 751,215 total 
nontrivial points for the 10,000,000 iteration). That is 
to say, near ten percents of nontrivial points need to be 
reconsidered when a third order SRA is taken instead of 
second order SRA. A natural conclusion is that one has 
to take the order of SRA as high as possible to ensure 
the validity of SRA. It seems that we are lucky enough 
today- as can be seen from Fig.7- the regions where Con- 
ditions B and C are violated seems to be disfavored by 
current observations. For the stringent condition A we 
can find the straight lines in Figs. 1,3 are due to the vio- 
lation of this condition. It is however noted that Condi- 
tions B and C are the weak conditions which work only 
to third order and nothing could be ensured under fourth 
or fifth order SRA. Would this be the real case, one has 
to solve Eq.l7 mode by mode instead. (A naive way out 
might be to ensure the validity of decreasing sufficiently 
the higher order fiow parameters.) It is worth mention- 
ing again that the fiow equations of Eq.9 is always exact 
under any order of truncation. In this sense the way 
of using flow equations like RefT^ is appropriate, while 
many other papers in the literature suffer from second or- 
der SRA0,[lI|l3ElEil. 

When a truncation to eighth 
order is considered instead we find similar results. For an 
iteration of 1000,000 points with truncation to the eighth 
order (Method Crosscheck) shown in Fig. 2 we get 1,972, 
289 and 5,216 points for models that violate Conditions 
A,B and C. A 10,000,000 iteration of Method Crosscheck 
is also tried and we get 17 insignificant points, 1,841 with 
I^AhI > 10*^ and 757,521 nontrivial points. The 19,719, 
2,804 and 51,626 points which violate Conditions A, B 
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and C are shown in corresponding window in Fig. 8. ^ 
This seems to put current inflationary flow to consider- 
able jeopardy. 

In summary, the WMAP result of a varying spectral 
index, if further stands, could be used as a discrimina- 
tor for inflationary models; most extant models would 
face a severe challenge. Using the inflationary flow equa- 
tions, we have studied in this paper the possibility of 
reconstructing inflation models with large running spec- 
tral indices, which is favored by the WMAP analysis. 
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Table 1 - Twelve sets of initial flow parameters that lead to large 
running in the index as shown in Fig. 4. 
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